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A.  Summary 


The  aim  of  this  project  was  to  develop  a  computational  modeling  capability  to  study  the 
theoretical  strengths  and  deformation  mechanisms  of  structural  ceramics  for  high  temperature 
applications.  The  approach  adopted  was  molecular  dynamics  simulation  based  on  appropriate 
interatomic  potential  models.  The  studies  focused  on  ZrC  which  has  excellent  high-temperature 
properties,  and  to  a  lesser  extent  on  SiC  and  Si02.  Strain-to-failure  simulations  were  performed 
to  observe  the  details  of  the  structural  instability  and  subsequent  microstructural  evolution. 
Besides  single  crystals,  deformations  in  the  amorphous  and  nanocrystalline  phases  were  also 
investigated.  In  parallel  elastic  stiffness  and  other  mechanical  properties  were  calculated  to 
correlate  with  and  interpret  the  simulation  results. 

The  three-year  grant  was  in  operation  effectively  over  four  years,  2000-2003,  as  a  result 
of  a  no-cost  extension  of  10  months.  The  major  thrusts  undertaken  during  each  year  were: 

(Yr.  1)  Developed  simulation  capabilities  for  strength  and  deformation  studies  involving 
hydrostatic  tension  and  compression  loading;  began  development  of  an  interatomic  potential  for 
ZrC.  (Yr.  2)  Continued  ZrC  potential  development  including  validation  tests;  began  construction 
of  a  model  of  silica  nanorod  for  tensile  failure  simulations.  (Yr.  3)  Performed  fracture  simulation 
studies  of  Si02  and  ZrC  using  interatomic  potentials  and  shear  deformation  calculations  using 
first-principles  density  functional  theory  codes;  continued  study  of  silica  nanorod.  (Yr.  4) 
Commenced  study  of  water-silica  reaction  using  quantum  chemical  techniques 

The  following  studies  are  considered  significant  accomplishments  of  the  project. 

1.  Development  and  validation  of  an  interatomic  potential  for  ZrC. 

2.  Theoretical  strength  of  beta-SiC  and  alpha-Si02 

3.  Uniaxial  compression  of  alpha-quartz 

4.  Fracture  toughness  of  crystalline  quartz 

5.  Shear  response  of  nanocrystalline  ZrC 

6.  First-principles  calculation  of  affine  shear  deformation  in  crystals 

7.  Shearability  gap  between  ceramics  and  metals 

8.  Silica  nanorod  model  for  deformation  and  reaction  studies 

9.  Kinetic  pathway  for  hydrolysis  of  strained  silica  nanorod 

Each  result  is  discussed  further  below,  followed  by  a  brief  section  on  conclusions  and 
recommendations. 


3 


B.  Research  Accomplishments 

1.  Development  and  Validation  of Interatomic  Potential  for  ZrC 

The  transition  metal  (refractory)  carbide,  ZrC,  is  known  for  its  hardness  (Vickers 
microhardness  of  25.5  GPa)  and  very  high  melting  point  (3700K),  and  high-temperature 
compressive  yield  strength.  Although  brittle  at  room  temperature,  it  undergoes  a  brittle-to-ductile 
transition  below  1300K.  In  order  to  model  its  mechanical  behavior  at  the  atomistic  level,  an 
interatomic  interaction  model  is  needed.  We  have  constructed  an  empirical  many-body  potential 
for  ZrC  with  a  functional  form  derived  from  the  formalism  of  Local  Combination  of  Atomic 
Orbitals  (LCAO),  a  so-called  tight-binding  method  in  electronic  structure  calculations.  The 
model  has  6  parameters  in  describing  the  Zr-C  bonding,  which  are  fitted  to  a  carefully  selected 
database  consisting  of  cohesive  energy  curves  for  three  cubic  structures,  elastic  constants  of 
NaCl-ZrC,  heat  of  formation,  and  most  importantly,  atomic  force  constants  accurately  determined 
by  first-principles  density  functional  theory  [1]. 

To  validate  our  model  we  systematically  calculated  the  following  properties  using  the 
potential  and  compared  them  with  either  more  rigorous  electronic-structure  results  or 
experimental  measurements  -  structure  and  energetics  of  hexagonal  WC  phase,  carbon  vacancy 
formation  and  migration  energies,  metal  vacancy  formation  and  migration  energies,  (001)  surface 
relaxation,  bond  shrinkage  near  carbon  vacancy,  Gruneisen  parameter  and  thermal  expansion 
coefficient,  vibrational  amplitude  and  thermal  strain,  and  the  melting  point.  None  of  these 
properties  were  used  in  fitting  the  parameters,  therefore  each  comparison  serves  as  an 
independent  test  of  the  robustness  of  the  model  description. 

We  find  our  potential  gives  a  better  description  of  the  various  properties  than  a  recent 
tight-binding  model.  The  reason,  we  believe,  is  that  the  hopping  integrals  in  the  latter  treatment 
do  not  have  any  angular  dependence.  Secondly  we  find  the  thermal  properties  of  ZrC  to  be  well 
described  by  our  model.  This  can  be  attributed  to  the  fact  that  the  model  gives  very  reasonable 
phonon  dispersion  curves  when  compared  against  experiments  and  theoretical  calculations. 
Fundamentally  this  goes  back  to  our  understanding  of  the  nature  of  bonding  in  ZrC,  namely,  a 
strong  bond-angle  dependence.  The  model  also  gives  temperature  variations  of  the  thermal  strain 
and  vibrational  amplitudes  in  good  agreement  with  measurements.  Direct  simulation  of  the 
melting  of  a  cluster  yields  3500K,  which  can  be  compared  to  the  experimental  value  of  3700K. 

2.  Theoretical  Strength  of  f  —  SiC  and  a  —  Si02 

Strength  calculations  have  been  performed  on  f  —  SiC  in  molecular  dynamics 
simulations  using  a  bond-order  potential,  a  purely  covalent  description  [2].  Cohesive  strengths 
have  been  determined  for  single  crystal,  amorphous,  and  nanocrystal  at  room  temperature.  The 
single  crystal  strength,  at  38  GPa,  represents  the  upper  limit.  The  amorphous  and  nanocrystal 


values,  being  close  to  each  other  (about  24  GPa),  suggest  a  certain  underlying  connection 
between  these  two  types  of  microstructures.  Similar  simulations  have  been  extended  to 
a  -  Si02  using  a  potential  description  consisting  of  short-range  and  electrostatic  interactions  [3]. 

Under  hydrostatic  tension  we  find  a  polymorphic  transition  from  a  to  P  phase.  Since  this  is  a 
well-known  transformation  when  quartz  is  heated,  we  have  demonstrated  the  equivalence 
between  thermal  and  mechanical  strains  [4].  Under  hydrostatic  compression  we  obtain  a 
transition  to  amorphous  at  20  GPa  [5]  which  is  close  to  the  experimental  value. 

3.  Uniaxial  Compression  of  Alpha-Quartz 

We  have  studied  alpha-quartz  under  uniaxial  compression  with  a  confining  pressure  of 
1.5  GPa.  This  mode  of  loading  is  chosen  to  correspond  with  experiments  in  geoscience  on  the 
longstanding  problem  of  hydrolytic  weakening  of  quartz  [6].  In  contrast  to  the  response  to 
hydrostatic  compression,  the  lattice  deformed  elastically  up  to  a  critical  stress  of  14  GPa,  and 
corresponding  strain  of  about  12%,  when  structural  instability  set  in.  Under  constant  pressure  the 
system  deformed  in  a  time-dependent  manner  on  the  scale  of  picoseconds,  showing  the  lattice 
becoming  disordered  in  a  heterogeneous  manner  of  nucleation  and  growth  of  a  planar  region.  To 
our  knowledge,  this  process  has  not  been  discussed  heretofore.  The  deformation  features 
revealed  by  our  simulations  have  correspondence  with  experimental  measurements  on  plastic 
deformation  of  quartz  single  crystals  [7].  The  data  showed  that  failure  occurred  by  sudden  rupture 
giving  rise  to  fractures  with  different  orientations  and  planar  structures.  We  propose  that  one  can 
distinguish  two  possible  mechanisms  for  planar  fracture  -  one  is  the  development  of  shear  stresses 
to  activate  certain  slip  systems  and  give  rise  to  material  failure,  and  the  other  is  the  onset  of  an 
elastic  instability  (under  compression)  to  nucleate  a  local  region  of  disorder  which  then 
propagates  through  the  lattice.  The  former  is  the  conventional  view  which  may  be  called  slip- 
assisted,  and  the  latter  is  the  one  we  are  presently  proposing  which  may  be  called  instability- 
assisted  [8].  By  adopting  the  instability-assisted  mechanism  of  failure,  it  is  possible  to  give 
plausible  explanations  of  several  reported  observations  [7]  that  are  not  consistent  with  the 
conventional  view  of  slip-based  mechanism. 

4.  Fracture  Simulations  of  Crystalline  Quartz 

As  part  of  our  investigation  of  strength  and  deformation,  molecular  dynamics  simulations 
of  crack  extension  of  under  mode  I  (uniaxial  tension)  loading  were  performed.  Two  systems 
were  investigated,  an  embedded  nanocrack  and  a  notch  at  a  free  surface  [3].  In  both  simulations 
the  crack  extension  occurs  at  a  characteristic  critical  stress,  lower  in  the  case  of  the  notch  as  is 
expected,  in  a  brittle  manner.  Contrasting  details  of  strain  localization  around  the  crack  tip  and 
bond  breaking  can  be  seen.  The  fracture  toughness  value  estimated  using  the  Griffith  criterion 
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compares  well  experiment,  while  the  value  obtained  using  the  Griffith  criterion  but  with  an 
unrelaxed  surface  energy  agrees  with  the  result  of  direct  simulation.  This  is  a  demonstration  that 
atomistic  simulation  and  experiment  probe  quite  different  time  scales  of  response. 

5.  Shear  Response  of  Nanocrystalline  ZrC 

We  have  deformed  several  small  nanocrystalline  samples  with  grain  size  varying  from 
2.4  nm  to  5.9  nm  by  applying  a  shear  stress  [3].  We  find  the  shear  strength  increases  with 
increasing  grain  size,  which  is  the  reverse  of  the  conventional  Hall-Petch  effect.  Previously  this 
behavior  has  been  reported  in  simulations  of  nanocrystalline  fee  metals;  our  results  indicate  that 
the  effect  is  quite  generic,  present  in  ceramics  as  well  as  in  metals. 

6.  First  Principles  Calculations  of  Affine  Shear  Deformation  in  Crystals 

Using  density  functional  theory  we  have  shown  that  aluminum  has  a  higher  ideal  shear 
strength  than  copper,  even  though  the  latter  has  a  higher  shear  modulus  [9].  This  result  is  a 
consequence  of  the  more  extended  range  of  deformation  before  softening  in  Al,  and  in  turn  it  can 
be  traced  to  the  directional  nature  of  the  chemical  binding.  By  comparing  the  ion  relaxations  and 
valence  charge  distributions  in  these  two  metals,  we  arrive  at  contrasting  descriptions  of  bonding 
characteristics  that  can  explain  their  relative  strength  and  deformation  behavior. 

7.  Shearability  Gap  between  Ceramics  and  Metals 

We  have  extended  our  first-principles  calculation  of  affine  shear  deformation  to  a  group 
of  fee  and  bcc  metals  and  covalent  and  ionic  ceramics  [10].  Focusing  on  the  maximum  shear 
strain  that  a  crystal  can  withstand  before  yielding,  we  find  spectral  gap  in  the  maximum  shear 
strain  between  ceramics  and  metals,  and  a  strong  correlation  of  this  shearability  with  the  degree 
of  directional  bonding.  Furthermore  we  show  that  incorporating  shearability  as  a  material 
parameter  one  can  capture  quantitatively  the  stress-strain  behavior  of  perfect  metallic  and  ceramic 
crystals,  thus  putting  our  fundamental  understanding  of  shear  deformation  on  the  same  footing 
with  the  universal  description  of  tensile  cohesion  in  solids.  Together  these  constitutive  relations 
allow  useful  materials  design  and  performance  criteria  to  be  formulated  where  shear  and  tensile 
deformation  modes  compete. 

8.  Model  of  Silica  Nanorod  for  Deformation  and  Reaction  Studies 

Current  interest  in  nanoscience  and  technology  has  stimulated  an  awareness  that 
mechanical  behavior  in  small  dimensions  can  have  significant  deviations  from  properties 
normally  associated  with  structures  at  the  two  extremes,  an  isolated  cluster  of  atoms  or  a 
macroscopic  specimen.  Unlike  physical  properties  which  deviate  from  continuum  behavior  only 
on  the  atomic  scale,  mechanical  strength  and  deformation  are  controlled  by  defects  in  the 
assembly  which  in  themselves  can  range  from  angstrom  to  nanometer  to  micometer  length  scales. 
We  have  introduced  a  nanorod  as  a  novel,  stand-alone  system  which  has  advantages  over  a  cluster 
in  that  meaningful  mechanical  behavior  can  be  extracted  and  scale  up  continuously  to  larger 
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dimensions.  The  rod  is  also  more  suitable  than  a  bulk  system  for  the  study  of  stress  and  chemical 
reaction  because  one  has  better  control  over  the  atomic  site  for  structural  defect  nucleation  or 
chemical  bond  reaction.  We  analyzed  mechanical  response  of  a  144-atom  nanorod  of  Si02  in 
uniaxial  tension  and  compression  and  compared  the  results  to  in-house  simulations  of  bulk  quartz 
[11].  This  work  is  a  prelude  to  the  study  of  the  effects  of  water  on  the  mechanical  strength  of 
silica,  which  is  a  fundamental  issue  in  understanding  the  dislocation  mechanism  of  hydrolytic 
weakening  [6]. 

9.  Kinetic  Pathway  for  Hydrolysis  of  Strained  Silica  Nanorod 

We  are  currently  investigating  a  particular  problem  of  environmental  effect  on  the 

mechanical  strength  of  a  ceramic  material,  that  of  hydrolytic  weakening  and  stress 

corrosion  of  silica.  We  are  interested  in  the  kinetics  of  hydrolyzing  the  strained  Si-0 

bond,  the  microscopic  mechanism,  the  reaction  pathways  and  rate.  That  the  presence  of 

water  has  a  significant  effect  on  strength  of  quartz,  particularly  at  high  temperature,  is 

well  known.  It  has  been  proposed  that  the  fundamental  mechanism  for  this  effect  is  the 

hydrolysis  of  a  siloxane  bond,  which  bridges  two  neighboring  Si04  tetrahedra,  to  form 

two  terminal  silanol  groups  [6].  The  silanol  groups  in  turn  will  facilitate  bond  rupture  and 

thus  lower  the  strength  of  silica.  Since  the  problem  involves  bond  breaking  and 

formation,  and  significant  effects  of  charge  transfer,  classical  potential  simulations  will 

not  be  adequate.  To  keep  the  system  to  a  size  that  permits  quantum  calculations,  we  are 

studying  the  nanorod  which  we  have  recently  developed.  Initially  we  are  applying  a 

semi-empirical  molecular  orbital  method  as  implemented  in  a  general-purpose  package 

MOP  AC.  Our  preliminary  results  on  sampling  the  reaction  pathways  indicate  a  three- 

step  mechanism  for  the  hydrolysis  --  adsorption  of  the  water  molecule,  water  reaction 

with  the  siloxane  bond  Si-O-Si  involveing  dissociation  and  proton  transfer,  and  formation 

of  silanol  group  SiOH  [12].  Knowing  the  pathways  we  are  able  to  determine  the 

activation  barrier  and  study  its  variation  with  strain  on  the  nanorod. 
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C.  Conclusions  and  Recommendations 

The  research  performed  in  this  project  has  demonstrated  the  feasibility  of  modeling  the 
intrinsic  mechanical  and  thermal  properties  of  ceramics.  The  results  obtained  contribute  to  a 
database  of  structure-property  correlation  which  provides  a  foundation  for  further  computational 
studies  of  ceramics  for  high  temperature  applications.  It  appears  that  modeling  and  simulation  at 
the  atomistic  level  have  reached  a  level  of  physical  robustness  that  one  can  begin  to  develop 
detailed  understanding  of  certain  cramics  under  conditions  of  practical  interest,  such  as  in  high- 
temperature  and  reactive  environments  in  synthesis  and  performance. 

Discussions  at  the  recent  Ultra  High  Temperature  Ceramics  Program  Review, 
Wintergreen,  VA,  November  5-7,  2003,  have  focused  on  a  number  of  issues  pertaining  to  the 
synethesis,  structural  evolution,  and  characterization  of  mechanical  and  thermal  properties  of 
diborides  and  carbides.  Specifically  there  is  a  need  for  a  more  fundamental  understanding  of 
strength  retention  at  high  temperature,  mechanisms  of  oxidation,  and  combustibility  at  high 
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oxygen  pressure.  It  is  our  belief  that  quantum  chemical  calculations  combined  with  atomistic 
simulations  of  the  kind  described  here  provides  a  capability  to  achieve  the  physical  understanding 
and  insights  that  can  help  to  interpret  existing  experimental  data,  as  well  as  guide  further  efforts 
in  synthesis  and  materials  development. 
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